Strengthening the Bell Theorem: conditions to falsify local realism in an experiment 
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The two-particle correlation obtained from the quantum 
state used in the Bell inequality is sinusoidal, but the stan- 
dard Bell inequality only uses two pairs of settings and not 
the whole sinusoidal curve. The highest to-date visibility of 
an explicit model reproducing sinusoidal fringes is 2/tt. We 
conjecture from a numerical approach presented in this paper 
that the highest possible visibility for a local hidden variable 
model reproducing the sinusoidal character of the quantum 
prediction for the two-particle Bell-type interference phenom- 
ena is ^. In addition, the approach can be applied directly 
to experimental data. 



It is a common wisdom in the quantum optical com- 
munity that the threshold visibility of the sinusoidal two- 
particle interference pattern beyond which the Bell in- 
equalities are violated is (for the case of perfect detec- 
tors) -T= (see, e. g. [Q). Most of the experiments exceed 
that limit (with the usual "fair sampling assumption") 
H. Some difhculties to reach this threshold were ob- 
served in the very early experiments m, as well as in 
some recent ones involving novel techniques. Thus far, 
in atomic interferometry EPR experiments H and for the 
phenomenon of entanglement swapping [^ , the resulting 
visibility is less than the magic 71%. 

It is also well known that the Clauser-Horne inequality 
and the CHSH inequality are not only necessary condi- 
tions for the existence of local realistic models but also 
sufficient M (in the case of the CHSH inequality, this re- 
quires some simplifying assumptions B). However, the 
sufficiency proofs used involve only two pairs of settings 
of the local macroscopic parameters (e.g. orientations 
of the polarizers) that define the measured local observ- 
ables. Thus, the constructions are valid for precisely 
those settings and nothing more, and there is no guar- 
antee that the models can be extended to more settings. 
Consequently one may ask what is the maximal visibil- 
ity for a model applicable to all possible settings of the 
measuring apparata, that returns sinusoidal two particle 
interference fringes. It is already known that for perfect 
detectors, this value cannot be higher than -7= or lower 
than 2/7r (this is the visibility of the recent ad hoc model 
by Larsson rah for earlier models returning visibilities of 
50% see e.g. [|). 

The knowledge of the maximal visibility of sinusoidal 
two-particle fringes in a Bell-type experiment that still 



can be fully modeled in a local realistic way, may help us 
to distinguish better between 'local' and 'nonlocal' den- 
sity matrices. For two two-state systems one can find 
precise conditions which have to be satisfied by density 
matrices describing the general state, pure or mixed, of 
the full system, that enable violation of the CHSH in- 
equalities p. States fulfilling such conditions are often 
called "nonlocal". However, since the CHSH inequality 
is necessary and sufficient only for two pairs of settings, 
it is not excluded a priori that some states that satisfy 
such inequalities for all possible sets of two pairs of lo- 
cal dichotomic observables, nevertheless give predictions 
that in their entirety cannot be modeled by local hidden 
variables. Such models must first of all reproduce the 
full continuous sinusoidal variation of two-particle inter- 
ference fringes, as well as the other predictions. 

It is clear that the full solution of the question would 
require a construction, or a proof of existence, of local 
hidden variable models which return sinusoidal fringes of 
the maximal possible visibility, that are applicable for all 
possible settings of the measuring apparata. Since this 
seems to be very difficult, we chose a numerical method of 
pointwise approximation at a finite number of settings at 
each side of the experiment. Due to exponential growth 
of the computation time when the number of settings in- 
creases, we managed to reach up to 9 settings on each 
side, i.e. up to 81 measurements points (which due to 
a certain symmetry, about which we will say more later, 
effectively can be transformed into 18 x 18 = 324 points). 
The exponential growth hinders any substantial increase 
in this number pC| . Such numerical models cannot give 
a definite answer concerning the critical visibility of si- 
nusoidal fringes, however our calculations enable us to 
put forward a strong conjecture that this value must be 
indeed 4= (see below). 

Experimentally our problem can be formulated in the 
following way: the two particle state produced by the 
source does not allow for single particle interference, and 
in the experiment less-than-perfect two-particle fringes 
are obtained, due to some fundamental limitations (like 
those present in the case of entanglement swapping, e.g. 
M) or due to imperfections of the devices. What is the 
critical two particle interference visibility beyond which 
the observed process falsifies local realism? We shall ask 
these questions assuming, for simplicity, perfect detec- 
tion efficiencies, which is possible theoretically, and ex- 
perimentally thus far amounts to the usual "fair sampling 
assumption" . 



Furthermore, the problem may be investigated with- 
out the use of the assumption that the observed fringes 
are of a sinusoidal nature even though two observed two- 
particle fringes in experiments with high photon counts 
follow almost exactly the sinusoidal curves [EJ . In exper- 
iments with lower count rates, still with relatively good 
level of confidence the recorded data have approximately 
the same character, and it is customary to fit them with 
sinusoidal curves. It is now a standard procedure to per- 
form the two-particle interference experiments by record- 
ing many points of the interference pattern, rather than 
stabilizing the devices at measurement settings appropri- 
ate for the best violation of some Bell inequality. Further, 
in some of the experiments, e.g. those involving optical 
fiber interferometers it is currently not possible to sta- 
bilize the phase differences and what is observed is just 
the interference pattern changing in time, and the visi- 
bility of the sinusoidal two-particle fringes is used as the 
critical parameter [[ll[ . 

Even though the numerical calculations presented here 
only reach 18 x 18 points, this is more than enough in 
comparison to the experimental data. The usual experi- 
mental scans rarely involve more than 20 points. Further, 
our algorithm can be applied directly to the measurement 
data, and in that way one can even avoid the standby hy- 
pothesis that the fringes follow a sinusoidal pattern. The 
algorithm can directly answer the question: are the data 
compatible with local realism or not? Since physics is 
an experimental science, the questions about Nature get 
their final answers solely in this way. 

Let us now go to a formal treatment of the problem, 
and our numerical solution of it. In a standard Bell- 
type experiment one has a source emitting two particles, 
each of which propagates towards one of two spatially 
separated measuring devices. The particles are described 
by the maximally entangled state, e.g.. 
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where |-|-)i is the state of the first particle with its spin 
directed along the vector jz of a certain frame of reference 
(— denotes the opposite direction), etc. We will now as- 
sume that the measuring devices are Stern-Gerlach ap- 
parata, measuring the observable n ■ cr, where n = a,b (a 
for the first observer, b for the second one), n is a unit 
vector representing direction at which observer n makes a 
measurement and cr is a vector the components of which 
are standard Pauli matrices. The family of observables 
n ■ (T covers all possible dichotomic observables for a spin 
2 system, endowed with a spectrum consisting of ±1. 

In each run of the experiment every observer obtains 
one of the two possible results of measurement, ±1. The 
probability of obtaining the result to = ±1 at the ob- 
server a, when measuring the projection of the spin of 
the incoming particle at the direction a, and the result 
I = ±1 at the observer b, when measuring the projection 
of spin of the incoming particle at the direction b is equal 
to 



while the probabilities of obtaining one of the results in 
the local stations reveal no dependence on the local pa- 
rameters, PQM{m\0') — \ and PQM{l]b) — ^- In a 
real experiment, however, one cannot expect that the 
observed probabilities will follow (g). Therefore, we will 
allow that the interference pattern is of a reduced visi- 
bility. In such a case (0) should be replaced by 



PQM{m, I; a,b) — j{l — mlVa ■ b), 



(3) 



where < V^ < 1 stands for the visibility. 

Now for the pointwise approximation. For V = 1 the 
quantum prediction for two-particle correlation function 
reads: 



EqAiia, b) = J^^^.i mlP{m, I; a, b) = -a ■ b, 



(4) 



and there is no single particle interference (i.e. the 
local results are absolutely random). If one assumes 
that the unit vectors which define the measured observ- 
ables are always coplanar the correlation function can 
be simplified to EQM{a,P) — —cos{a — /3) (with the 
obvious definition of a and (3). Let us assume that in 
the experiment the observer at the side a chooses be- 
tween, say, N settings of the local apparatus, denoted 
here by a^ with i = 1,2, 3,... TV, and the other ob- 
server, at side 6, chooses between M settings bj, with 
j = 1,2,..., M. The quantum correlation function for 
V — 1 has at these settings the following fixed numerical 
values: EqM{o.i,bj) = —ai-bj. Thus, we have a cer- 
tain matrix of quantum predictions. We can denote this 
matrix by Q, with Qij — EQMi<ii,bj). For V < I the 
correlation function reduces to VQij. 

The Bell theorem is a statement on the impossibility 
of modeling certain quantum predictions by local and 
realistic theories. To make the matters simpler let us 
use the local hidden variable (LHV) formalism. Within 
such a formalism the correlation function must have the 
following structure 



ELHviai,bj)^ / dXp{X)A{ai,X)B{bj,X), 
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where for dichotomic measurements A{ai,X) = ±1 and 
B{bj,X) = ±1, and they represent the values of local 
measurements predetermined by the LHV's, denoted by 
A, for the specified local settings. This expression is an 
average over a certain LHV distribution p(A) of certain 
factorizable (rank 1) matrices, namely those with ele- 
ments given by M[X)ij — A{ai,X)B{bj, X). The sym- 
bol A may hide very many parameters. However, since 
the only possible values of A{ai, A) and B{bj, A) are ±1 
there are only 2^ different sequences of the values of 
{A{ai, X), A{a2, X), ..., A{aN , X)), and only 2'^' different 
sequences of {B{bi,X),B{b2,X), ...,B{bM, A)) and conse- 
quently they form only 2^+*^ matrices M(A). Therefore 
the structure of LHV models of ELnvidi, bj) reduces to 



discrete probabilistic models involving the average of all 
the 2^+*^ matrices A/(A). In other words, the LHV's 
can be replaced, without any loss of generality, by a cer- 
tain pair of variables k and I that have integer values 
respectively from 1 to 2^ and from 1 to 2*^. To each 
k we ascribe one possible sequence of the possible val- 
ues of A(ai,X), denoted from now on by A{ai,k), sim- 
ilarly we replace B{bj,X) by B{bj,l). With this nota- 
tion the possible LHV models of the correlation function 
EhRvio-i^bj) acquire the following simple form 



ELHv{ai,bj) ^ ^^pkiA{ai,k)B{bjJ), 



(6) 



fc=i 1=1 



with, of course, the probabilities satisfying pki > and 

J2Pki = 1- 

The special case that we study here enables us to sim- 
plify the description further. To satisfy the additional 
requirement that the LHV model returns the quantum 
prediction of equal probability of the results at the local 
observation stations, that P{1; a) = P{m; b) = i, one can 
use the following observation. For each k, there must ex- 
ist a fc' ^ fc with the property that A{ai, k') = —A{ai, fc), 
and similarly for each I, there must exist an /' ^ / for 
which B{bj,l') = -B{bjJ). Then A{a^,k)B{bj,l) = 
A{ai,k')B{bj,l'), and thus they give exactly the same 
matrix of LHV predictions. By assuming pki — Pk'i' the 
property of total randomness of local results will always 
be reproduced by the LHV models, and the generality 
will not be reduced since the contributions of pti and 
Pk'i' to (0) cannot be distinguished. Of course in the ac- 
tual computer calculations of the correlation function we 
take only one representative of the two pairs, reducing 
in this way the number of probabilities and matrices of 
LHV predictions in (g) by a factor of two. 

Another reduction by a factor of four is given by the 
fact that in the coplanar case, the choice of the settings 
may be limited on each side to ranges not greater than 
TT (i.e. (j) < ai < (j) + -K,^ < f3j < ij) + it). This is due 
to the simple observation that a model of the type (||) 
once established for such settings can be easily extended 
to settings a^ — aj 4-7r, /3- = j3j+tt by putting A(a^, I) — 
— A(ai,/) and B{(3j,l) = —B{(3j,l) (nevertheless, some 
scans with wider ranges were performed). 

The conditions for LHV's to reproduce the quantum 
prediction with a final visibility V can be simplified to the 
problem of maximizing a parameter V for which exists a 
set of 2^+^^ probabihties, pki, such that 



^ ^ Pfez A(a„ fc)B(bj , = ^g^. 



(7) 



fc=i 1=1 



This is a typical linear optimization problem in which we 
have more unknowns than conditions and for which many 
good algorithms exist. Therefore the simplest method to 
solve it is to use the standard method of linear program- 
ming. The core of the algorithm is a procedure which 



finds the maximum of a linear function within the given 
constraints. In our case the constraints are the N x M 
equations (|7|), and the condition that J^Pki = 1- Our 
unknowns are all pki and V (all nonnegative) . We can 
treat them as points in 2^^^ + 1 dimensional space. 
The constraints given by (0) define some subset of this 
space. On this subset we use the trivial linear function 
f{pkh T^) = T^ as our goal function (for clarity, our func- 
tion depends only on the variable V) and we search for 
its maximum. 

In order to apply this method to experimental results, 
replace Qij in (^) by the measured values E°^^ , and per- 
form the same task. If the critical V ||l^ returned by the 
program is less that 1, the data cannot be reproduced 
by any LHV model. Note that one even does not have 
to know what the settings are(!). This type of approach 
may be useful especially when one is not able to stabi- 
lize the interferometers at settings which are required for 
some Bell inequality, but there are data at many other 
settings. 

For the case of coplanar settings we have checked lots 
of "interesting" combinations of the settings of the ap- 
parata at each side (e.g. N x N problems , 2 < A^ < 9, 
with evenly spaced settings, etc. [ffsl). The main inter- 
esting generic feature of our results is the following one. 
Whenever among the coplanar settings a^ and /3j there 
is a subset ai^ , ai^ and pj^ , Pj^ such that for these set- 
tings the CHSH inequality (equivalently CH inequality) is 
maximally violated by the ideal quantum prediction, our 
optimum, the maximal visibility reproduced by LHV's, 
is 1/ = -^. In all other cases, namely for those settings 
without any such a subset, we have obtained maximal 
visibilities describable by LHV's (usually, slightly) higher 
than -7= . However with increasing number of more or less 

evenly spaced settings this difference decreases. For gen- 
eral measurement directions (i.e. those including non- 
coplanar settings) all numerical scans follow the same 
pattern as for the co-planar case. Additionally, we have 
checked 200000 randomly chosen sets of 5 x 5 settings, 
and never a visibility lower than -y= has been returned. 
All this strongly suggests that a LHV model returning 
correlation function of the form -y=a ■ b in its entirety 
indeed exists. 

Furthermore, in the coplanar case, using carefully cho- 
sen equally spaced settings more structure is introduced 
into the model, in the form of two additional symmetries 
of the matrix Q (in addition to the ones described above) : 
ordinary matrix-transpose symmetry, and constant diag- 
onals. This reduces the rate of the exponential growth 
of the calculation to a point where the problem is com- 
putable on a standard PC in reasonable time, even for 
13 X 13 (extendible to 26 x 26) settings. Even in this case 
the returned visibilities exhibit the behavior discussed 
above and always satisfy V > -y=. 

Finally let us present application of our method to raw 
experimental data. In a recent Bell-type experiment We- 
infurtcr and Michler have obtained the following matrix 



of results |14| 










' -0.894 


-0.061 


0.761 




-0.851 


0.343 


0.765 




-0.625 


0.688 


0.516 




-0.251 


0.860 


0.103 


E^p = 


0.226 
0.530 


0.921 
0.651 


-0.389 
-0.648 




0.855 


0.323 


-0.832 




0.852 


-0.092 


-0.843 




0.785 


-0.539 


-0.638 




0.397 


-0.795 


-0.253 



(8) 



Our program gives the verdict that the values of all en- 
tries to the matrix of results have to be reduced by the 
factor of 0.796 to be describable by local hidden variables. 
In the recent long-distance EPR-Bell experiment the 
following set of values of the correlation function was 
obtained [1151: 



^cxp _ 



0.960 


-0.102 


0.903 


-0.375 


0.733 


-0.660 


0.479 


-0.809 


0.191 


-0.903 


-0.120 


-0.923 


-0.429 


-0.807 


-0.666 


-0.656 


-0.842 


-0.395 


-0.951 


-0.152 


-0.953 


0.171 



(9) 



This matrix has to be reduced by the factor of 0.7366 to 
have a local realistic description. 

To conclude, the performed calculations enable us to 
put forward the following conjecture, which is the main 
result of this letter: sinusoidal two-particle fringes of vis- 
ibility up to -7= are describable by local realistic theories. 
At this stage we are not able to give an analytic proof of 
the above. However, for finite sets of measurement points 
the results of data analysis with the use of our program 
fully concur with this hypothesis. This implies, e.g. that 
one needs a re-run of the entanglement swapping exper- 
iment, in order to show that this phenomenon can lead 
to observable violations of local realism. 
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